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Summary. — The outflows from gamma ray bursts, active galactic nuclei and rel- 
ativist ic jets in general interact with the surrounding media through collisionless 
shocks. With three dimensional relativist ic particle-in-cell simulations we investi- 
gate such shocks. The results from these experiments show that small-scale mag- 
netic filaments with strengths of up to percents of equipartition are generated and 
that electrons are accelerated to power law distributions N(j) oc j~ p in the vicinity 
of the filaments through a new acceleration mechanism. The acceleration is locally 
confined, instantaneous and differs from recursive acceleration processes such as 
Fermi acceleration. We find that the proposed acceleration mechanism competes 
with thermalization and becomes important at high Lorentz factors. 
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By the way, which PACS is it, the 00.00? GOK. 



1. — Introduction 

Observations of astrophysical objects with relativist ic plasma outflows generally show 
frequency spectra dominated by power-law segments. These objects count gamma ray 
bursts (GRBs) and active galactic nuclei. The general interpretations of such observations 
are that 1) electrons are accelerated to power-law distributions N(j) oc j~ p in the 
collisionless shock interface between the relativistic eject a and the surrounding media 
and 2) the accelerated electrons radiate through synchrotron radiation in the downstream 
region of the shock. The magnetic field strength required to account for the observed 
radiation spectra is typically of the order of percents of equipartition. A complete physical 
model for relativistic collisionless shocks should account for both the origin of such strong 
magnetic fields and for acceleration of electrons to power-law distributions. 

Regarding the origin of the magnetic field, particle-in-cell (PIC) simulations have 
shown that the Weibel two-stream instability can generate a small scale, but strong, 
magnetic field transverse to the jet flow (cb — 0.05 where e# describes the fraction of 
total injected kinetic energy that is converted to magnetic energy) [IJ 03 E3 E3 C3 • Other 
PIC simulations have found that the ambient large scale magnetic field neccesary to 
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Fig. 1. - Left panel: The Weibel two-stream instability forms current channels. The electrons 
(top panel) undergo the instability faster the the heavier ions (bottom panel). Right panel: 
A schematic view of an ion current channel surrounded by an electric and a magnetic field. 
Electrons in the vicinity of the current channels are subject to a Lorentz force with both an 
electric and magnetic component, working together to accelerate the electrons along the ion flow. 
Crossing the center of the channel, the process reverses, leading to an oscillating movement along 
the channel. 



quench this instability has to be significantly stronger than the /aG field typically found 
in the interstellar medium (ISM) j^l E] . 

Regarding the non-thermal electron acceleration, Monte Carlo test-particle simula- 
tions have shown that Fermi acceleration can provide electron power-law distributions 
under some assumptions about the shock and magnetic field 8 . This mechanism has, 
however, not been conclusively demonstrated to occur in ab initio particle simulations. 
PIC simulations have shown that both thermal and non-thermal particle acceleration 
take place in the shock transition region [2J EI EH • 

In this paper, we discuss the particle dynamics in collisionless shocks based on results 
from a three dimensional relativistic PIC code. We propose a new particle acceleration 
mechanism that is itself a consequence of the Weibel two-stream instability 0. 

2. — Numerical Model and Results 

We have performed computer experiments using a charged particle-in-cell code. The 
code works from first principles by solving Maxwell's equations for the electromagnetic 
fields and the Lorentz force equation of motion for the particles. The fields are defined 
on a numerical grid and the particles are defined in a continuous phase space. We use a 
computational box with 125 x 125 x 2000 grid points and a plasma consisting of 8 x 10 8 
ions and electrons. The ion rest-frame plasma frequency is co P i = 0.075, rendering the 
box 150 ion skin depths long. The electron rest-frame plasma frequency is co pe = 0.3 in 
order to resolve the microphysics of the electrons. Hence, the ion-to-electron mass ratio 
is rrii/m e = 16. 

In the experiments we study the microphysics of two colliding plasma populations. 
The experiments are carried out in the rest frame of one of the populations (downstream, 
e.g., a jet). In this frame, a less dense population (upstream, e.g. the ISM) is continuously 
injected at the leftmost boundary z = 0. We have performed two runs with relativistic 
velocities corresponding to Lorentz factors T = 3 and T = 15. Initially, the jet plasma 
is denser by a factor of three and both plasma populations are unmagnetized. The 
boundaries in the direction transverse to the jet flow (x,y) are periodic. At the z- 
boundaries, electromagnetic waves are absorbed and we allow particles to escape in order 
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to avoid unphysical feedback. The experiment last until t max = 340 , which is 
sufficient for the continuously injected upstream ISM-plasma to travel 2.3 times the 
length of the box. 

The results from the experiments show how the Weibel two-stream instability forms 
current filaments in the region where the two plasma populations interpenetrate and 
how these filaments merge into increasingly larger patterns further downstream (Fig. 
left panel) (see also [I] E2])- The induced magnetic field grows to e# — 1%. We 
identify two mechanism that accelerate electrons. The first is thermalization caused by 
random deflections by the generated small scale electromagnetic field. This mechanism 
is dominant for the T = 3 run. In the T = 15 run a second acceleration mechanism 
becomes important. Here we briefly discuss this mechanism (for details please see 0): 

The generated ion current channels are Debye shielded by the heated electrons. At 
distances less than the Debye length, the current channels are surrounded by transverse 
electric fields that accelerate the electrons toward the current channels. The acceleration 
happens in an approximately electrostatic field and is a simple consequence of potential 
energy being converted into kinetic energy. Therefore, the electrons are decelerated again 
when leaving the current channel and reach their maximal velocities at the centers of the 
current channels (Fig. Q}. It has been shown that a spatial Fourier decomposition of the 
transverse ion filaments in the box exhibits power-law behavior Hence, the number 

of accelerated electrons is expected to reflect this power-law behavior. This relation is 
confirmed in Fig. The left panel shows a scatter plot of electron momentum (try) 
plotted against the strength of the electric current sampled at each electrons position. 
There is a power-law correlation between the electron energy and the local current density. 
Some deviations from the power-law is found: 'Cold' trapped thermal electrons (indicated 
with the ellipse) exist inside the ion current channel and count towards lowering the 
average four velocity at high J ion . Also, heated electrons are found outside the current 
channels which increase the average four velocity at low Jion- The electron distribution 
function can be seen in the right panel of Fig. It show a thermal distribution with a 
power-law component at high energies. The slope of the power-law segment corresponds 
to a spectral index of p=2.7. 

3. — Conclusions 

We have investigated magnetic field generation and particle acceleration in collision- 
less shocks using three dimensional charged particle-in-cell experiments. The results are 
applicable to interactions between relativistic outflows and the interstellar medium. Such 
relativistic outflows occur in GRBs and in jets from compact objects. 

We observe how the Weibel two-stream instability is excited and find that the non- 
linear stage of the instability penetrates at least 150 ion skin depths into the shock ramp. 
Extrapolating to the real electron-ion mass ratio, this is more than 6000 electron skin 
depths. The instability efficiently creates a strong, small-scale magnetic field (cb — 1%). 

Regarding particle acceleration, the results show that thermalization dominates the 
electron distribution in the simulation run where a jet is expanding with T = 3. At 
T = 15, a new acceleration mechanism becomes important in the high energy tail of 
the spectrum. The acceleration is non-thermal and the resulting electron distribution 
function has a power-law segment. The slope of this segment corresponds to p = 2.7. 

We emphasize that the generation of magnetic field and the non-thermal acceleration 
of electrons are two highly interconnected processes in collisionless shocks. The acceler- 
ation is closely tied to the spatial distribution of electric current filaments. Further large 
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Fig. 2. - Left panel: A scatter plot of the local ion current density Ji on versus the four velocity 
of the electrons in a region downstream of the shock. Overplotted is a line (thin) showing the 
average four velocity as a function of Ji on , and a line (thick) showing a straight line fit. Right 
panel: The normalized electron particle distribution function downstream of the shock. The 
dot-dashed line is a power law fit to the non-thermal high energy tail, while the dashed curve 
is a Lorentz-boosted thermal electron population. The horizontal line at vj = 15 indicates the 
plasma injection velocity. 



scale simulations are need to determine how this slope depends on different parameters. 

When the injected ISM ions has traveled through the simulation box they are not 
merged with the jet ions to a single population. This indicates that larger simulations 
are needed to cover the whole shock ramp. 

* * * 
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